Inflating with the QCD Axion 
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We show that the QCD axion can drive inflation via a series of tunneling events. For axion models 
with a softly broken Zn symmetry, the axion potential has a series of N local minima and may 
be modeled by a tilted cosine. Chain inflation results along this tilted cosine: the field tunnels 
from an initial minimum near the top of the potential through a series of ever lower minima to the 
bottom. This results in sufficient inflation and reheating. QCD axions, potentially detectable in 
current searches, may thus simultaneously solve problems in particle physics and provide inflation. 



In 1981, Guth lj proposed an inflationary phase of 
the early universe to solve the horizon, flatness, and 
monopole problems of the standard cosmology. During 
inflation, the Friedmann equation 
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is dominated on the right hand side by a (nearly con- 
stant) false vacuum energy term p ~ p vac ~ constant. 
The scale factor of the Universe expands superluminally, 
a ~ t p with p > 1. Here H = a/a is the Hubble pa- 
rameter. With sufficient inflation, roughly 60 e-folds, the 
cosmological shortcomings are resolved. 

Standard inflationary models require the invention of 
a new field whose potential drives the superluminal ex- 
pansion; there is no direct evidence that the associated 
particle exists. In this paper, we demonstrate that such a 
new field is unnecessary: it is possible for the QCD axion 
a to drive inflation. The QCD axion has been proposed 
|3, 0] as a solution for the strong CP problem in the 
theory of strong interactions. It is advantageous to use a 
single particle to solve several problems. The mass scales 
of the axion are much lower than those of standard infla- 
tionary models, and the axion may be found in ongoing 
experiments, especially axion searches. 

While the axion is a priori a Goldstone boson of the 
spontaneously broken Peccei-Quinn symmetry U(1)pq, 
QCD instanton effects induce an axion potential with 
residual Zn symmetry. Our model includes an additional 
explicit soft-breaking term, which tilts the instanton in- 
duced potential. While the complete form of the axion 
potential is dependent on non-perturbative effects, it is 
well modeled as shown in Figure 1 by 
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The first term models the periodic instanton potential 
as a cosine with N degenerate vacua, or N bumps. The 
width of each bump is given by the Peccei-Quinn scale 
/„ = V /N ~ (10 9 - 10 12 ) GeV, and the height of each 
bump Vo — to 2 ./, 2 ~ QCD scale. The second term in (J5J 
is the tilting effect of the soft-breaking term. 
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FIG. f : Potential energy density of the QCD axion field 9 = 
a/ f a as a function of 0. The soft-breaking potential is a tilted 
cosine as in Eq. (J2j. Here we have taken N = 20 and n — Vo. 



Inflation starts with the axion field located in a mini- 
mum at the top of the tilted cosine potential. The uni- 
verse tunnels to the next minimum in the cosine, then 
on down through all the minima until it reaches the bot- 
tom. The universe inflates a fraction of an e-fold while 
it is stuck in each of these minima. Sufficient inflation 
results for N ~ few hundred. The general framework of 
a sequential chain of tunneling fields was considered pre- 
viously in the Chain Inflation model proposed by two of 
us 0- 

We consider the invisible axion model of Zhitnitskii 
and Dine, Fischler, and Srednicki (DFSZ) [1 0. The 
axion is identified as the phase of a complex SU(2) x U (1) 
singlet scalar a below the PQ symmetry breaking scale 

where a — -^{ v + p) ex P(*^)- The periodicity of 
the axion field is a = a+2nv = a+2nNfa, where f a is the 
axion decay constant. Defining 9 = a/ f a , we see that 8 is 
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2kN periodic. Below the QCD scale A QC d ~ 220MeV, 
QCD instantons produce a potential with N degenerate 
minima at 9 = 2ixn where n = 0, 1, 2, . . . , N — 1. The 
effective action is 



£ ff = ^d^ad^a 



(m 2 Jt)g(a/f a ), 



(3) 



where g(x) is a periodic function of period 2-7T, Taylor 
expanded as g(x) — g(0) — (l/2)a; 2 +•••. The axion 
mass is m a ~ 2N f^m^ / v , where z = m u /m,d — 0.56. 
Then m a f a ~ m^f^ where the pion decay constant 
/,r = 93MeV and the pion mass = 135MeV. Here N 
refers to the unbroken Zjq subgroup of U(1)pq, and cor- 
responds to the number and representations of fermions 
that carry color charge as well as PQ charge (and thus 
contribute to the QCD anomaly) (Jj. Since we need 
N ~ 200, this introduces additional heavy fermions be- 
yond the usual quarks and leptons. The form of g{x) 
depends on non-perturbative effects, and hence is not 
fully specified. We will take it to be g(x) — cos(x). This 
captures the main features of the periodic instanton po- 
tential, and will be sufficient for our purposes. 

We will take the Upq(1) symmetry to be softly broken 
• Following , we add a soft breaking term of the form 
£soft = A* 3 "" + h.c. where /i is a complex parameter and a 
is given above. Below the PQ scale, this adds a term to 
the axion potential of the form r\ cos(a/i>+7) where r\ and 
7 are real parameters. The combined potential is then 
given in Eq. J2J. Note that the phase shift 7 misaligns 
the QCD and soft breaking minima. Without loss of 
generality, we may restrict 7 to lie in the range —tt/N < 
7 < tt/N, Away from the bottom of the potential, the 
tilt can be treated in the linear regime, so that the total 
potential is of the form 
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The energy difference between minima is roughly e ~ 
2ttt]/N. In this linear regime, the requirement that e be 
less than the barrier height becomes e < Vq. Unless this 
criterion is satisfied, the barrier becomes irrelevant and 
the field simply rolls down the hill. 

Sufficient Inflation and Reheating. A successful infla- 
tionary model has two requirements: sufficient inflation 
and reheating. Here we have a series of a large number 
(TV) of tunneling events as the universe transitions from 
an initial high vacuum energy down to zero. In order 
to have sufficient inflation, the universe expands by at 
least 60 e-folds once all the tunneling events have taken 
place (by the time the field travels all the way down the 
cosine): 



Xtot > 60, 



(5) 



where Xtot is the total number of e-folds. In order for the 
universe to reheat, the number of e-folds attained dur- 
ing one tunneling event must be small (less than 1/3 of 



an e-fold), as shown below. The failure of old inflation, 
known as the "graceful exit" problem, is circumvented 
because the bubbles of vacuum are able to percolate at 
each step down the potential since the phase transition is 
fairly rapid. Chain Inflation's basic mechanism of mulit- 
pile tunneling events 0] works both in the context of the 
stringy landscape or here with the QCD axion. 

In the zero-temperature limit, the nucleation rate T 
per unit volume for producing bubbles of true vacuum in 
the sea offalse vacuum through quantum tunneling has 
the form 



T(t) = Ae 



(6) 



where Se is the Euclidean action and where A is a de- 
tcrminantal factor which is generally the energy scale e 
of the phase transition}^ . 

Guth and Weinberg have shown that the probability of 
a point remaining in a false deSitter vacuum is approxi- 
mately 



4-7T 

p(t)~ex.p(-—pHt), 
where the dimensionless quantity (3 is defined by 
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Writing Eq. Q as p(t) ~ exp(— t/r), we estimate the life- 
time of the field in the metastable vacuum as roughly [23j 
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The number of e-foldings for the tunneling event is 
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The authors of EH and [12 calculated that a critical 
value of 



9/47T 
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is required to achieve percolation and thermalization. In 
terms of e-foldings, this is 



X < Xc 



1/3. 



(12) 



As long as this is satisfied, the phase transition at each 
stage takes place quickly enough so that 'graceful exit' is 
achieved. Bubbles of true vacuum nucleate throughout 
the universe at once, and are able to percolate. 

Tunneling Rate. In the thin wall limit, the tunneling 
rate is given by Eq. JSJ . As shown by H, @ ; we nee d to 
calculate 



Si = J V2U+(a) da, 



(13) 
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integrated from one minimum to the next, where the 
symmetric portion of the potential is 



U+{9) = V (l-cos< 



Then 



Si = y^m/a / Vi-cosfl de = 8/„ 
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The Euclidean action is 



5 = 
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For the parameters of the DFSZ axion, So ^ 1 and tun- 
neling is suppressed in the thin wall limit (n.b. the thin 
wall limit almost never applies to any realistic tunneling 
event for any potential as tunneling is suppressed |l3|'). 
To have reasonably fast tunneling with % < 1/3 from one 
minimum to the next, we must be outside the thin wall 
limit. There is the additional constraint e < Vq in order 
for tunneling to take place, as opposed to mere rolling 
down the potential. Thus, obtaining the right amount of 
inflation requires e/Vb ~ 1/2. 

The Neutron Electric Dipole Moment. We must ensure 
that the soft-breaking term in Eq. (J2J does not destroy 
the strong CP solution, i.e., that the minimum of the 
potential in Eq. @ is not shifted away from zero by more 
than is allowed by the electric dipole moment (EDM) of 
the neutron [ljl 



Atf| EDM <6xlO 
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(17) 



To find the minima of the potential, we solve V'(a) = 0, 
or 

V Nsin(Na/v) + r]sin(a/v + j) = 0. (18) 

To leading order in small r\, the minima are located at 
a n = 2nnf a — ?7^y sin(2^- + 7) for integer n where the 
potential is V(a n ) — — 7/cos(2^ +7). The energy differ- 
ence between two adjacent minima is 
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The difference in field value between minima is 



Sa = 2nf a +i] 
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We have to impose the EDM bounds at the bottom 
of the potential, at n — 0, since this is presumably the 
endpoint of tunneling (corresponding to the current uni- 
verse). For large N, we find that the shift from = is 
given by 
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In the last equality, we have used that I7I < ir/N to 
estimate that a typical arbitrary value of 7 ~ it/ (27V). 

During most of the route down the potential, away 
from the bottom, the tilt can be approximated as being 
linear, as in Eq. J3J. In the linear regime, e = 2irr]/N. 
Taking, e.g., n ~ N/4, and using Eq. i|21|). we find 



4NV M\ 



EDM' 



(22) 



Combining this with the bound on the neutron EDM, 



we find that en, 



: 2 x 10" 9 iVA 



QCD: 



,1/4 



01 e Uncar < 



5(7V/200) 1 / 4 MeV. To get a sensible reheat temperature 
Tr > lOMeV after inflation 24] where Tr ~ eij ncar , we 
need a large number of new heavy fermions. 

At the bottom of the potential, e(n = 0) ~ 2ir 2 ri/N 2 , 



ig Eq. HQ). e bottom - 47rl/ A6»| 



EDM' 



Combining 



this with the bound on the neutron EDM, we find that 

bottom < 2Me V- 



The two conditions in the two regimes 
e < kK% cd M\ 



I EDM' 



(23) 



where k = 4-7T at the bottom of the potential and k — 
AN in the linear regime. Substituting this equation into 
Eq. (|16fl . we see that the Euclidean action can be written 
as 



r> x hi ' [ A 

. m„ 



1 



k 3 (A8 



EDM 



(24) 



The tunneling rate at the last stage is extremely sup- 
pressed for parameters allowed by the constraints on the 
neutron EDM, 0|edm < 6 x 10~ 10 . To obtain a reason- 
able tunneling rate, we need to get away from the thin 
wall limit (as discussed previously); i.e. the value of e 
must be closer to Vq ~ (lOOAfeV) 4 and hence must be 
larger than allowed by Eq. I|23fl ■ In order for Chain Infla- 
tion with the axion to succeed, we must reconsider some 
assumptions we have made. 

With large enough AT, the tunneling rate can be per- 
fectly reasonable as long as one stays away from the ab- 
solute bottom of the potential [see Eq. (|2^|l with n = AN 
in the linear regime] . Indeed, the reheating of the uni- 
verse can take place in the linear regime. As the field 
goes farther down the potential, the vacuum energy gets 
smaller and smaller, and fewer e-folds result. Hence, ra- 
diation that is produced during reheating that takes place 
•near, but not at, the bottom of the potential, is not in- 
flated away by the last few episodes of tunneling to the 
bottom. We have not yet investigated details of the par- 
ticles produced during reheating; axions may provide the 
dark matter. 

The only real problem with the model is the last 
tunneling event, near the bottom of the potential. In 
Eq. Ij23(l , with n = Air near the bottom, we see that there 
is no dependence on A^ in this regime. The allowed en- 
ergy difference is simply so small that no tunneling takes 
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place. The universe would still be situated in this false 
vacuum now. 

Circumventing the Neutron EDM constraint: The con- 
straint on the neutron EDM prevents the universe from 
tunneling in the last stage to zero energy. We may spec- 
ulate about a number of ways to resolve this problem. 
It is possible for the tunneling to stop at an energy of 
^bottom < 10~ 3 eV and thereby account for the dark en- 
ergy. Alternatively, the soft breaking term and the QCD 
axion may conspire to set 7 = in, avoiding the EDM 
constraint altogether. Third, one might consider a time 
dependent tilt or a different function for the tilt. Alter- 
natively, coupling to other fields may allow the universe 
to tunnel or roll to the bottom via a different direction 
in a two dimensional potential in field space. Or if the 
Zjf symmetry breaks once the field is near the bottom 
of the potential, then the field could quickly roll to zero 
energy. Further afield, axion models different from the 
DFSZ axion (see, e.g., the review of might work 
better as inflaton candidates. Briefly, there are models 
with several axions |l7lll8l . li9| with more desirable prop- 
erties 25j. Axions abound in string theory. Some may 
solve the strong CP problem 0, [3]. We have listed 
a few approaches here to circumvent the neutron EDM 
constraints on tunneling to the bottom. 

On the Value of N. The value of N need not equal 
the number of fermions that carry color and CP charge. 
The number of fermions may be far less, depending on 
the relevant group representations. Defining to be 
the generators for the representation r, one can write 
Tr(T ( " if J = \t T b a P. The axion model has an exact 

Zn symmetry with N = (2n/Te){J2 r ,i QnU) where Q r i 
are the PQ charges of the fermions and Tq, the period of 
O, need not be 2tt (see e.g. the review of Sikivie 21]). 

Domain Walls. One might worry about the delete- 
rious effect of domain walls which appear when differ- 
ent horizon sized portions of the universe fall into dif- 
ferent minima at the QCD scale. As shown by Sikivie 
, the energy difference between vacua leads to pressure 
that causes the domain walls to move in the direction 
of eliminating the higher energy vacua in favor of the 
lower energy ones. The domain walls disappear within a 
Hubble time. One might worry that the universe is too 
quickly driven to zero energy, without any inflation, but 
this does not happen. By the Kibble mechanism, there 
are at most two or three domains in any horizon volume, 
with different values of (a). Typically a horizon volume 
will be driven to a field value half way down the poten- 
tial, at ~ N/2, and will subsequently inflate sufficiently. 
In some part of the universe, a horizon volume will con- 
tain only domains with values of (a) ~ 0, near the top of 
the potential. Those regions that start the highest up the 
potential will inflate the most. Hence, if one plunks down 
our observable universe in a random patch of the total 
universe after inflation ends, then one is likely to find a 



region that started inflating near the top of the poten- 
tial and hence inflated sufficiently. Domain walls even 
have the positive effect of driving a causal patch prior to 
inflation to become more uniform by shoving away the 
nonuniformities . 

Conclusion. We have investigated using the QCD ax- 
ion potential to inflate. We use the cosine shape of the 
axion, with N minima, due to a residual Zjv symmetry, 
together with a tilt produced by a small soft breaking 
of the Peccei Quinn symmetry. We studied the DFSZ 
axion. Chain inflation results along this tilted cosine, 
with a series of tunneling events as the field tunnels its 
way from a minimum near the top of the potential to 
ever lower minima. Sufficient inflation as well as reheat- 
ing result. Tunneling in the last stage is suppressed due 
to constraints on the neutron Electric Dipole Moment 
and must be further considered. We have listed a few 
attempts to work around these constraints. 

In this paper, we have restricted discussion to axions 
which can solve the strong CP problem. Obviously, if 
we forego any contact with real QCD, then the allowed 
ranges for parameters becomes much larger. For ex- 
ample, the constraint from the neutron EDM vanishes. 
Then the ranges of potential width, barrier height, and 
energy difference between vacua are completely opened 
up. A tilted cosine may arise due to (non-QCD) "axions" 
in many other contexts, such as string theory, and would 
easily provide an inflaton candidate. 

We thank A. Guth, H. Haber, G. Kane, P. Nilles, P. 
Sikivie and E. Weinberg for useful conversations. We 
thank the Aspen Center for Physics as well as the Michi- 
gan Center for Theoretical Physics for hospitality and 
support while this work was completed. This work was 
supported in part by the US Department of Energy under 
grant DE-FG02-95ER40899. 



* Electronic address: ktfreese@umich.edu 
t Electronic address: jimliu@umich.edu 

* Electronic address: dspolyar@physics.ucsc.edu 
[1] A. H. Guth, Phys. Rev. D23, 347 (1981). 

[2] K. Freese and D. Spolyar, Chain Inflation in the Land- 
scape: Bubble Bubble Toil and Trouble, JCAP 0507, 
007 (2005). 

[3] S. Weinberg, Phys. Rev. Lett. 40, 223 (1978). 

[4] F. Wilczek, Phys. Rev. Lett. 40, 279 (1978). 

[5] A. P. Zhitnitskii, Yad. Fiz. 31, 497 (1980) [Sov. J. Nucl. 

Phys. 31, 260 (1980)]. 
[6] M. Dine, W. Fischler, and M. Srednicki, Phys. Lett. 

B104, 199 (1981). 
[7] P. Sikivie Phys. Rev. Lett. 48, 1156 (1982). 
[8] C. Callan and S. Coleman, Phys. Rev. D16, 1762 (1977). 
[9] S. Coleman, Phys. Rev. D15, 2929 (1977); M. B. 

Voloshin, I. Yu. Kobzarev, and L. B. Okun, Yad. Fiz 

20, 1229 (1974) [Sov. J. Nucl. Phys. 20, 644 (1975)]. 



5 



[10] S. Coleman and F. de Luccia, Phys. Rev. D21, 3305 
(1980). 

[11] A. H. Guth and E. Weinberg, Nucl. Phys. B212, 321 
(1983). 

[12] M.S. Turner, E.J. Weinberg, and L.M. Widrow, Phys. 

Rev. D46, 2384 (1992). 
[13] E.J. Weinberg, private communication. 
[14] P. G. Harris et al, Phys. Rev. Lett. 82, 904 (1999). 
[15] T. Banks and M. Dine, Nucl. Phys. B505, 445 (1997). 
[16] T. Banks and M. Dine, Phys. Rev. D68, 075011 (2003). 
[17] S. M. Barr, X. C. Gao and D. Reiss, Phys. Rev. D26, 

2176 (1982). 

[18] K. Choi and J. E. Kim, Phys. Lett. B132, 342 (1983). 
[19] J. E. Kim, H. P. Nilles and M. Peloso, JCAP 0501 005 
(2005). 



[20] J. E. Kim, Phys. Rept. 150, 1 (1987). 

[21] P. Sikivie, UF-TP-83-6 Based on lectures given at 21st 

Schladming Winter School, Schladming, Austria, Feb 26- 

Mar 6, 1982. 

[22] We note that we do not need to include gravitational 
effects |10J as they would only be relevant for bubbles 
comparable to the horizon size, whereas the bubbles con- 
sidered in this paper are much smaller. 

[23] There exists a distribution around this typical value. 

[24] This is the minimal requirement to obtain ordinary ele- 
ment abundances from nucleosynthesis. Then baryogen- 
esis must take place later, as considered by [Te|1 . 

[25] In superstring Es x E$ models, two axions are present. 



